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The dynamics of an extremely diluted neural network with high 

order synapses acting as corrections to the Hopfield model is inves- 

^ ' tigated. As in the fully connected case, the high order terms may 

strongly improve the storage capacity of the system. The dynamics 
displays a very rich behavior, and in particular a new chaotic phase 
emerges depending on the weight of the high order connections e, 
the noise level T and the network load defined as the rate between 
the number of stored patterns and the mean connectivity per neuron 
a = P/C. 
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1 INTRODUCTION 

Neural networks have been subject of intense research in Statistical Mechan- 
ics in the last decade. However, besides the equilibrium properties that are 
reasonable well understood in the framework of spin glass theory, dynamical 
properties are hard to treat analytically and still requires further work. In 
the limit of low loading or focusing only in the first few time steps, several 
treatments [ij yielded important results, the challenge remaining for satu- 
rated networks 0. Particularly, for extremely and asymmetrically diluted 
networks the dynamics can be exactly solved 0. It is worth to stress here 
that both dilution and asymmetry are biological realistic ingredients absent 
in the fully connected Hopfield model. 

In order to introduce asymmetry in the connections, the synapses Jij and 
Jji are cut with probability 1 — C /N independently of each other [Q, where 
C is the mean connectivity of each neuron and the extreme dilution limit is 
obtained when C -C In A^ {N is the size of the network). 

Measuring the network load by a = P/C, where P is the number of stored 
patterns, Derrida, Gardner and Zippelius (hereafter DGZ) found that there 
is a critical value a^ = 2/tt above which the system cannot retrieve the stored 
information. Since the storage capacity is here measured by the amount of 
bits that can be stored (PN) per synapse (CN), in the diluted version the 
Hopfield model is more efficient than the fully connected case (C = N) in 
which Oe^ 0.138 @]. 

Synapses connecting more than two neurons can be considered to both im- 
prove the storage capacity of the network [Q-|TT| and to mimic real synapses 



existing in nervous systems (see P] and references therein). For instance, 
axon-axon-dendrite connections can be described as third order synapses 
(even more intricate connections, involving more than two axons, may also 
exist in the brain) [^]. Moreover, when some pairwise connections are close 
enough, they may interact somehow, and that can also be considered as high 
(even) order synapses. Diluted networks with high order connections were 
studied in refs. |0, |, 0] . Kanter and Tamarit et al considered a mono- 
mial of order k in the overlap and found a discontinuous transition at ac{k) 
for k > 2. Also, in the (a,T) phase diagram the retrieval phase shrinks as 
k increases: a^k) —>■ a.s k —>■ oo although the retrieval is perfect (m — >■ 1). 
Wang and Ross treated instead a polynomial, controlling the relative co- 
efficients (weights) and found, besides a retrieval phase, regions where the 



system can be either periodic or chaotic, depending on the noise or the degree 
of parallehsm in the updating. 

In this paper we study the dynamical effects of introducing a correction 
term to the Hopfield model [|Tl| in the extreme dilution limit of both, two 
and fourth order synapses. The difference between our model and those 
studied by Kanter et a/ 0, Tamarit et al Q and Wang et al lays in the 
nature of high order connections we consider (see details below). The paper 
is organized as follows: section |^ reviews the model and in sections ^ and ^ 
analytical and numerical results are presented for the dynamics. In section 
1^ we summarize and present our conclusions. 

2 THE TRUNCATED MODEL 

In a binary neural network, as for example the Hopfield model |T2|, |13[, each 
neuron is modeled by an Ising variable Si that can take the values {+1,-1} 
representing the active and passive state, respectively. Possible states of 
the network are given by A^-dimensional vectors S G {+1,-1}^ and the 
embedded memories are associated to P of these states, denoted by ^^ (/i = 
1,...,P). 

The fully connected Truncated model can be regarded as the Hopfield 
model plus correction terms. Since the synapses are symmetric, the dynam- 
ics is ruled by a Lyapunov function (called Hamiltonian by analogy with 
magnetic systems) and the long time behavior of the network can be inferred 
from a thermodynamical analysis of the system. In this case the Hamiltonian 



is given by |T0|, |TT 

i=\ fii<...<fj.i 

Ml Ml<M2 /ll<...</ii\f 

where 2M is the highest order of interaction and the overlap m^ between the 
state of the network S and the pattern ^^^ is given by 



m 



-1 "" 
1=1 



E^^S, ■ (2) 



In what follows we consider only the first correction to the Hopfield term 
(M = 2 and 821 = 5u + £^2£) • Eq- (H) can then be rewritten as 

'^^ = — - 2^ JijSiSj + - 2^ JijkiSiSjSkSi . (3) 

The learning rule for the second order Jij and fourth order Jijki follows di- 
rectly from eq.(|I]) and will be detailed in the next section [|13|. 



The equilibrium statistical mechanics ^ of this model was studied in 
refs. [|lT|, showing an extremely rich behavior that depends on the value of 
e. The main features of the system at T = (for non zero temperatures the 
behavior is analogous) are summarized by the phase diagram in fig.|l|: for large 
e (~ 0.5) the overlap m decreases continuously with a, going monotonically 
from 1 down to zero at a'^{e)] as e decreases (~ 0.36), the overlap presents 
a local minimum, before finally going to zero at a^{6). If e < ^c — 0.3587, 
the minimum yields a gap separating two retrieval regions. Not too close to 
Ec the gap A has a discontinuous border at a'^ and a continuous one at a~ . 
The points where m continuously approaches zero, af{e), are 




^c{e)={^±\- . (4) 



As e — i> 0, A = a^(£) — o^ci^) goes to infinity as e^^. Thus, when the Trun- 
cated model recovers the Hopfield model (e — *> 0), the location of the second 
retrieval region in the a-axis goes to infinity and a^(0) ~ 0.138. Numerical 



simulations [|Tl|] at T = show the above transitions and the existence of the 
gap. The basins of attraction seem to be large and a-independent; the mean 
convergence time increases with a and, for non small values of a, does not 
depend on the initial overlap. 

There is also an optimal value of e that strongly improves the storage 



capacity of the model by canceling the noise term for each value of a [^ 



^opt-Y^ . (5) 

As we will see in the next two sections, there is an analogue of this optimal 
value for the diluted network. A complete discussion of these and other 



results for the fully connected Truncated model can be found in refs.|]TT 



3 THE DILUTED NETWORK 

In this section we study the dynamics of an extremely diluted and asymmetric 
version of the network described in section ^ in two different cases: the 
initial state has a macroscopic overlap with a) only one pattern and b) two 
correlated patterns (out of P — 2 uncorrelated ones) . Due to the asymmetry 
in the connections, a Lyapunov function can no longer be defined for the 
network and hence we are constrained to study the time evolution of the 
system ruled by the Heat Bath dynamics, given by 

f +1 with probability (1 + exp[-2pohi{t)])-'^ 
Siit + At) = <^ (6) 

[ -1 with probability (1 + exp[+2 (3ohi{t)])-^ 

where the parameter jSo = T^^ (called the inverse of the temperature) mea- 
sures the noise level of the net and hi{t) is the local field acting on the neuron 
i at time t: 

hi = Y^ JijSj - sY^ J-jj,iSjSkSi . (7) 

j j,k,l 

Based on eq.(|I|), the couplings are: 



■Jijki - o ("'^i'^' + ■^^J''i + "^'^i*') (^) 



1 

3 

Jijki = Cijki 22 ^i^j^k^i (10) 

and Cij and Cijki are random variables distributed according to the proba- 
bilities p and p: 

p(a,) = §'5(C., - 1) + (l - ^) 5(a,) 



(11) 



C f c 

p{Ctjki) = j;p^{Cijki ~ ^) + (^ ~ jv3 ) ^C^vki) 



where C <^ log A^. The asymmetry is introduced through the independence 
of Cij and Cji (the same holds for Cijki). 



In the fully connected case some of the self-connections are not zero, par- 
ticularly, Juki and Jijkk- These self-interactions create correlations between 
the states of neuron i at different times, that is, {Si(t)Si(t')) ^ 0, what appar- 
ently prevents from using the DGZ prescription to obtain the time evolution 
of the diluted network. Actually, it can be proved that the field generated 
by the self-interactions vanishes in the thermodynamical limit, implying null 
correlations. 

Assuming that the updating is parallel and the initial state is correlated 
with only one of the embedded memories, we are interested in obtaining a 
recurrent equation, 

m(t + l) = /(m(t)) , 

for the overlap mif) between the state of the network and this memory, 

^W = ^E(C^^^W) ' (12) 

i 

where (...) denotes both thermal and configuration averages. Using standard 
techniques and considering non biased patterns, after taking the limit of 
C — > cxo, the equation ruling the parallel evolution of the system reads: 



Vy tanh /3 m(t) - V2ai/(1 - em^{t)) , (13) 

where Ca = P — 1, /3 = C jT^ and T>y is the gaussian measure: 

Vy=%e-^' . (14) 



We also study the case where the initial state has a macroscopic overlap 
only with the two first memories, which have a fixed overlap between them, 

^ = \iie.e. , (15) 

i 

the remaining P — 2 being uncorrelated. The quantities of interest are then 
the overlaps with these two memories: 



^i(t) = ^E(e-^.W) 



^2(t) = ^E(e-^.W) 



(16) 



and the time evolution equations for them are ruled by: 

M(t + 1) = (I + k) fVytanhp 



Mil- e- - y\^ 1 - ^ -^ - 



(17) 



m 



m{t + l) = {1 - k) Vy tanh (3 

where: 

M = nil + 1^2 

m = mi — 7712 



(M^-m^)\ ^( eiM^ + m') 
1 + e -yV2^[l 



(18) 



In the next section the above equations are numerically solved and the 
results with their interpretation are presented. 

4 RESULTS 



In the T = limit, fim) becomes 



/(m) = erf 



m{t) 



2a{l-em^{t))_ 



(19) 



Fig.|^ displays this map for two different values of e showing that for £ < 1 
it is continuous while it is discontinuous otherwise. 

In figj^ we show the fixed points of the equation m{t + 1) = f{m{t)) as 
a function of a for several values of e < 1. For e < 0.266 the system has 
only a fixed point that decreases continuously down to at Oc = 2/11. Above 
this critical value, there is only the m = solution. For 0.266 < e < 1 the 
transition is discontinuous and ac is greater than 2/7r, increasing with e. In 
this case there are three different regimes: for a < 2/tt there is only the 
retrieval solution; for 2/tt < a < ac this solution coexist with the m = one 
and for a > ac there is only the paramagnetic solution. In the intermediate 
regime, an unstable solution appears that separates the basins of attractions 
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of the two stable solutions and the retrieval basin decreases with a as can be 
seen in fig.^. Also, for a given value of a, the closer e is to 1 the better is the 
retrieval. The retrieval region below the line Oc, where exists a nonzero fixed 
point, is named Ri, while above a^ the paramagnetic phase is named P . 

At e = 1 there always is the solution m = 1 and the unstable solution 
found for a > Ijv: drops to zero as we approach e = \. Hence, the basin 
of attraction of the retrieval solution increases while the basin of the param- 
agnetic solution shrinks to zero when e — i> 1. This means that the network 
capacity diverges and the retrieval is perfect for all nonzero initial states! 
This behavior is the analogue of the Eopt found in the fully connected system 
where e = (1 + a)~^ while here we have 

_ 1 1 _ 

^opt - YTpJn - 1 + aC/N ~ ^ ' ^^°^ 

since C <^ \nN . 

For values of e > 1 a rich behavior emerges and the system presents a 
novel kind of retrieval. The overall behavior of the system is presented in 
fig.^, showing several regimes. For a fixed e and low values of a the system 
is in a retrieval phase (R2) where it oscillates between a state with overlap 
m with the pattern and a state with —m overlap. 

As a increases the system enters in a new phase. It starts oscillating 
between two distinct values of m, say mi and m^ (|?tii| 7^ \'m2\) while higher 
periods appear with increasing values of a leading to a chaotic regime. The 
route to chaos followed by the system is not exactly period doubling since 
it suffers split bifurcations |T^: before doubling the period of the attractor. 



the system doubles the number of stable attractors. For instance, a period 
2 cycle splits in two period 2 cycles before becoming a period 4 one. The 
transition from the twofold attractor to the double period one occurs when 
the system enters a superstable orbit (those that contain the critical points). 
This behavior is very similar to that found in the cubic logistic map f{m) = 
am{l — m^) where only one split bifurcation is allowed in the periodic region 
(a system with two critical points can only have two stable attractors) |I^] . 
A representative behavior of the system is shown in fig.^. The basins of each 
attractor for a = 0.6 and e = 2 were also examined: for an initial value 
of m around the domains are extremely mixed. This behavior is better 
understood if one plots the function iterated several times (e.g., 16) as a 



function of m (fig. H): the function /^^^^ is a set of discontinuous plateau at 
some values of m. 

To decide whether a given temporal sequence is chaotic or not we calcu- 
lated the Lyapunov exponent through its definition 

1 " 
A= lim-^ln|/'K)| (21) 

"- i=l 

where n is the length of the sequence. In the aperiodic regime one has A > 
indicating that the system is chaotic while if its behavior is periodic, A < 0. 
It is worth noting that even in the periodic and in the chaotic attractors, the 
system always passes very close to the memorized pattern, as can be observed 
in figj^ for e = 2. We have verified numerically that for any point in this 
phase (C), the attractors have at least one point with overlap m near 1, 
hence it is possible to interpret this cyclic or chaotic phase as an alternative 
retrieval phase, in which the system does not recognize a memory by reaching 
a fixed point but by wandering around it. 

For large enough a the system enters a paramagnetic phase (the only 
fixed point of eq. ([l9D ) as can be observed in the phase diagram fig.^ Notice 
also the similarities between the T = phase diagram for the diluted and the 
fully connected system: both have a retrieval region around the optimal value 
of 6 for all values of a due to the canceling of the noise term. The difference is 
that in the fully connected model Eopt = Sopt{ci) while in the diluted one it is 
a constant. Also, in the fully connected case there is not a chaotic phase (not 
even a periodic one) because there the connections are symmetric allowing 
the introduction of a Lyapunov function (the Hamiltonian eq.(|l])). We have 
also determined that, as £ — ;> 1 



^) 



a,^\l-e\-' . (22) 

Figure ^ displays the T versus a phase diagram for e < 1 where there 
are only fixed point attractors: one with m = (P) and another with m ^ 
(Ri). In particular, for values of e where the transition is continuous the line 
Tc is independent of e and satisfies 

Vy sech^ (3cyV2^ = 1 • (23) 

-oo '- -' 



If e > 1, the map eq.([T3|) presents a solution with m = ±e^/^ at 

/5-f .n§±j . (24, 



It means that for T = T* = (/3*)^^ the system presents the above solution 
for all values of a. In the limit e — > 1, T* = and m = ±1, as is shown 
in fig.^ and, as e — > oo, T* — > 1. In fig.|^ the T versus a phase diagram for 
e = 2 is shown. Note that, besides the three phases appearing at T = 0, we 
now found also an Ri phase even for e > 1. For increasing values of a the 
m = region dominates except for points around T*. In fig.|^ the behavior 
of m for e = 2 and a = 1 is shown as a function of T with the corresponding 
Lyapunov exponent A. Note the backward period bifurcation leading to a 
retrieval phase for high temperatures: the system leaves the chaotic regime 
and become periodic as the temperature increases until a certain temperature 
where it has a fixed point and is able to retrieve the stored information. 

Finally we have analyzed the case in which only two memories have macro- 
scopic overlap k, being the dynamics ruled by eq.(|l8D, starting always from 
the initial condition mi(0) = 1 and m2{0) = n. The attractors of these equa- 
tions were numerically studied and four different phases were distinguished: 
a retrieval one (Ri) which correspond to a fixed point with \mi\ > \m2\ > 0, 
another retrieval phase (R2) in which the system oscillates in a cycle two 
between the states with |mj| = cte, \mi\ > \m2\, a mixed phase (M) with 
mi = 1712 7^ 0, a paramagnetic phase (P) with mi = 1712 = and a cyclic or 
chaotic one (C). 

For a = 0, when one considers macroscopic overlap with only one pattern, 
the DGZ equation is recovered. But, in the situation where the system has 
macroscopic overlap with two memories we do not recover the DGZ equations, 
that is, there is still an e dependence on the equations. Surprisingly, this also 
holds for K = 0: the system may present a cyclic or chaotic regime, absent if 



one looks for a retrieval solution. In fig. 10 we present the phase diagram T 
versus e for k = 0.2 and a = 0. For T = the system present only two phases: 
the retrieval one for \e\ < k,~^ and the cyclic and chaotic one for \e\ > k,~^, 
while the mixed phase appears for a given temperature that depends on the 
value of K. It is important to stress that the boundary between the retrieval 
phase and the cyclic or chaotic one for positive e shows the first appearance 
of cyclic orbits: inside the C phase there are islands of retrieval. Also, in 
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the boundary between the mixed phase and the paramagnetic one, there is a 
point £*, above which the transition is continuous and discontinuous below. 
It can be shown that the transition hne in the continuous case is given by 

T^ = l + K. 

In the case a 7^ and T = 0, the a versus e phase diagram is shown in 
fig.|ll]. This phase diagram is qualitatively similar to the one where there is 
overlap with only one memory (fig.§). There is also an value of e where the 
capacity diverges, but here it is given by 

£opi = l-fi:2 . (25) 

Besides that, the boundary between the paramagnetic phase and the mixed 
one is 

a,= -{l + Kf . (26) 

TT 

In the line a = of R2 region on the above diagram, there is always the 
retrieval solution Ri. 



5 CONCLUSIONS 

We presented an exact solution for the dynamics of a model for neural net- 
works with high order interactions that shows a very rich behavior depending 
on the parameters a,e and T. In analogy with the related fully connected 
model, there is an optimal value of e for which the model always retrieves 
the embedded information (although the process may not be perfect). It is 
important to stress that for some values of e the retrieval region is found 
at high values of T: amazingly, as the temperature is decreased, the system 
passes from an disordered phase (P) to an ordered one (R) and, after passing 
through a region where it presents cyclic or chaotic orbits, it reenters in the 
paramagnetic region. In this case, the presence of thermal noise may improve 
the retrieval abilities of the system. 

The periodic or chaotic behavior is only present for e > 1. Apparently, 
the reason for this is that below e = 1 the Hopfield second order term surpass 
the fourth order one while for e > 1 it is the reverse. In the later case, due to 
the asymmetry on the synapses (the presence of both yU and v in the learning 
rule), cycles are introduced. This is analogous to some models where the task 



is to retrieve temporal sequences [0 and the learning rule explicitly stores 
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the transitions between states. The absence of this kind of transition in the 
fully connected model is due to the possibility of symmetrization of all the 
connections, and hence allowing the introduction of an energy function for 
the system. Also, if e < 1 the map eq.(|l^) is invertible and continuous, not 
allowing chaos for one dimensional systems |l^ . 



Due to the very complex nature of the system behavior, it would also be 
interesting to study the thermodynamics of the symmetric diluted case 0], 
trying to see the interplay between the fully connected case and the highly 
diluted one. Also, if chaos can be controlled as proposed by Ott, Grebogi and 



Yorke |T8[, the phase boundaries of the above phase diagrams may change, 
increasing the retrieval region. 

The presence of the mixed (M) phase is an indication that the study of 
the generalization abilities of the model may reveal some interesting features. 

At last, a remark on the possible role of chaos on the problem of per- 



ception and recognition is worth. It was proposed |T^ that the property 



of large groups of neurons changing abruptly their activity pattern due to 
small inputs is responsible for the flexibility of real nervous systems to create 
new patterns when interacting with the outside world, flexibility that may 
underlie the learning and creative processes. 



Acknowledgments: We acknowledge R.M.C. de Almeida for several 
interesting discussions and for encouraging this work. Also, we are in debt 
with J. CM. Mombach, F.B. Rizatto, D. Stariolo and C Tsallis for useful dis- 
cussions. F.T.A. acknowledges the kind hospitality of the IF-UFRGS during 
his stay at Porto Alegre. Work partially supported by brazilian agencies Con- 
selho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq), Finan- 
ciadora de Estudos e Projetos (FINEP) and Fundagao de Amparo a Pesquisa 
do Estado do Rio Grande do Sul (FAPERGS). 



12 



References 

[1] COOLEN A.C.C. AND RuiJGROK Th.W. 1988 Phys. Rev. A38 4253 
Kepler T.B. and Abbott L.F. 1988 J. Phys. France 49 1657 
Gardner E., Derrida B. and Mottishaw P. 1987 J. Phys. France 
48 741 

[2] Patrick A.E. and Zagrebnov V. A. 1991 J. Phys. A24 3413 

[3] Derrida B., Gardner E. and Zippelius A. 1987 Europhys. Lett. 
4 167 

[4] Amit D.J., Gutfreund H. and Sompolinsky H. 1987 Ann. Phys. 
NY 173 30 

[5] Peretto p. and Niez J.J. 1986 Biol. Cybern. 54 53 

[6] Gardner E. 1987 J. Phys. A 20 3453 

Abbott L.F. and Arian Y. 1987 Phys. Rev. A36 5091 
Horn D. and Usher M. 1988 J. Phys. France 49 389 

[7] Kanter L 1988 Phys. Rev. A38 5972 

[8] Tamarit F.A., Stariolo D.A. and Curado E.M.F. 1991 Phys. 
Rev. A43 7083 

[9] Wang L. and Ross J. 1991 Phys. Rev. A44 R2259 

[10] DE Almeida R.M.C. and Iglesias 1990 Phys. Lett. A 146 239 
Arenzon J.J., DE Almeida R.M.C. and Iglesias J.R. 1992 J. 
Stat. Phys. 69 385 

[11] Arenzon J.J., de Almeida R.M.C., Iglesias J.R., Penna T.J. P. 
AND DE Oliveira P.M.C. 1993 Physica A 197 1 
Arenzon J.J. and R.M.C. de Almeida 1993 Phys. Rev. E48 4060 

[12] HOPFIELD J.J. 1982 Proc. Natl. Acad. Sci. USA 79 2554 

[13] Hertz J., Krogh A. and Palmer R.G. 1991 Introduction to the 
Theory of Neural Computation (Addison- Wesley Publishing Company) 



13 



[14] Testa J. and Held G.A. 1983 Phys. Rev. A28 3085 

[15] Skarda C.A. and Freeman W.J. 1987 Behavioral and Brain Sci- 
ences 10 161 

[16] For a review see KuHN R. AND VAN Hemmen J.L. in "Physics of 
Neural Networks" (E. Domany, J.L. van Hemmen and K. Schulten eds.), 
Springer 1990. 

[17] Canning A. and Gardner E. 1988 J. Phys. A21 3275 

[18] Ott E., Grebogi C. and Yorke J. 1990 Phys. Rev. Lett. 64 1196 

[19] Ott E. 1993 Chaos in Dynamical Systems (Cambridge University 
Press) 



14 



Figure Captions 



Figure 1: Phase diagram a versus e at T = of the fully connected Truncated 
model. The lines (x!: are second order while a' is discontinuous. 



Figure 2: Right hand side of eq.(]T^) for a = 0.5 and two values of e: 0.5 
(dashed) and 2 (solid). 



Figure 3: Fixed points m versus a for several values of e. The solid line is 
for stable solutions while the dashed is for unstable ones. For a > 2/n the 
paramagnetic solution is always stable. 
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Figure 4: Phase diagram a versus e at T = for the diluted Truncated 
model. The transition is continuous for e < 0.266. The optimal value of e is 
1 where Oc — > cxd as |1 — £:|~^. For e > 1 the system can present a periodic 
or chaotic phase (C) and a retrieval one (R). 

Figure 5: Plot of m versus a at T = for e = 2. For high values of e there 
is only the m = solution while for lower values of e the system behavior is 
complex. 



Figure 6: Map eq.(0) iterated 16 times showing the existence of several 
plateau for e = 2 and a = 0.6. 



Figure 7: Phase diagram for some values of e < 1 showing the existence of 
only two phases: a retrieval one (T < Tc) and m = one (T > Tc). For 
e < 0.266 the line T^ is the same as of the case e = 0. 



Figure 8: Phase diagram for e = 2. The retrieval {R) and the periodic or 
chaotic phases (C) are surrounded by the m = one (P). Note that for 
T* ~ 0.8 the system has the solution m = ±l/v^ for all values of a. 



Figure 9: The overlap m versus T for e = 2 and a = 2 and the correspond- 
ing Lyapunov exponent. The system has only a retrieval regime at high 
temperatures. 

Figure 10: Phase diagram for a = and k = 0.2 in the case where there is 
macroscopic overlap with two memories (see text). 

Figure 11: Phase diagram for T = and k = 0.2 in the case where there is 
macroscopic overlap with two memories (see text). Below the R2 phase, at 
a = 0, the system is in an Ri phase, and that below M, in a C phase. This 
compatibilizes this diagram with the former one (fig.p!0|). 
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